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The cellular recognition of viruses evokes the secretion of type-I interferons (IFNs) that induce an 
antiviral protective state. By live-cell imaging, we show that key steps of virus-induced signal 
transduction, IFN-p expression, and induction of IFN-stimulated genes (ISGs) are stochastic events 
in individual cells. The heterogeneity in IFN production is of cellular — and not viral — origin, and 
temporal unpredictability of IFN-p expression is largely due to cell-intrinsic noise generated both 
upstream and downstream of the activation of nuclear factor-KB and IFN regulatory factor 
transcription factors. Subsequent ISG induction occurs as a stochastic all-or-nothing switch, where 
the responding cells are protected against virus replication. Mathematical modelling and 
experimental validation show that reliable antiviral protection in the face of multi-layered cellular 
stochasticity is achieved by paracrine response amplification. Achieving coherent responses 
through intercellular communication is likely to be a more widely used strategy by mammalian cells 
to cope with pervasive stochasticity in signalling and gene expression. 
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Introduction 

The type-I interferon (IFN) system provides a powerful 
defence against viral infections (Kunzi and Pitha, 2003; 
Sadler and Williams, 2008; Takeuchi and Akira, 2009). 
Pathogen recognition receptors, such as RIG-I, stimulate an 
intracellular signalling cascade that leads to the activation of 
nuclear factor kB (NF-kB) and IFN regulatory factors (IRFs) 3 
and 7 (Brennan and Bowie, 2010) . Subsequent induction of the 
IFN- (3 gene and other type-I IFNs is a hallmark of the early 
response to infection (Theofilopoulos et al, 2005). Upon 
secretion and binding of type-I IFNs to their specific receptor, 
the Jak/STAT signalling pathway is activated to reprogram 
gene expression (Brierley and Fish, 2005) . By activation of a 
wide set of these genes, IFNs act directly by establishing an 
antiviral state and indirectly through recruiting innate and 
adaptive immune cells. 

On the one hand, IFN production must be tightly regulated 
to avoid harmful inflammation and autoimmune disease 
(Trinchieri, 2010). On the other hand, pathogenic viruses 
inhibit the expression of IFNs or IFN-stimulated genes (ISGs) 
by diverse mechanisms, indicating that maintaining the 
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efficiency of the IFN system in the face of these viral challenges 
might be an overriding objective of its evolution (Versteeg and 
Garcia-Sastre, 2010) . Therefore, the finding that only a fraction 
of cells in a virus-infected cell population expresses IFNs has 
been surprising (Zawatzky et al, 1985; Hu et al, 2007). Bio- 
chemical studies have suggested that stochastic IFN induction 
results from host cell-intrinsic causes, such as a mechanism of 
IFN gene induction (Hu et al, 2007; Apostolou and Thanos, 
2008) or cellular heterogeneity in expression of the viral sensor 
RIG-I (Hu et al, 2011). Alternatively, heterocellularity of IFN 
expression could be caused by the infecting virus (Chen et al, 
2010; Killip et al, 2011). It is conceivable that several factors 
might shape the dynamics of IFN induction, depending on host 
cell type and virus. Live imaging provides a suitable tool to 
study dynamics and variability at single-cell resolution (Spiller 
et al, 2010), but, such an analysis has not yet been carried out 
for the IFN system. 

To understand the functional consequences of cell-to-cell 
variability in IFN induction, the cellular response to IFNs, 
particularly the expression of ISGs, must also be defined. Two 
recent quantitative analyses of IFN-stimulated signalling have 
modelled the dynamics at the cell-population level (Maiwald 
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et al, 2010) and considered effects of cell-to-cell variability 
(Levin et al, 2011). However, the expression of ISGs and the 
resulting cell fate have not yet been characterized at single-cell 
level. 

Here, we study in living cells and at single-cell resolution 
both type-I IFN induction and the cellular response to secreted 
IFN. Using fluorescent reporters based on bacterial artificial 
chromosomes (BACs) and chimeric transcription factors 
(Supplementary Figure SI), we imaged successive key steps 
of IFN induction and response in a prototypical model system, 
the infection of murine cells in culture with a single-stranded 
RNA virus. To link our single-cell data to antiviral protection at 
the cell-population level, we developed a mathematical model 
based on these data and tested its predictions experimentally. 

Our results show that cell-to-cell heterogeneity is a 
pervasive feature of the IFN system. This heterogeneity 
manifests itself not only in the virus-induced expression of 
IFN but also in the IFN-induced protective response. It is to a 
large part due to cell-intrinsic stochasticity occurring in three 
key steps: virus-induced signalling, IFN gene induction, and 
expression of ISGs. In contrast to gene-expression noise 
described in bacteria and yeast that causes protein levels to 
fluctuate around a mean value in individual cells (Raj and van 
Oudenaarden, 2008), the stochasticity found here in the IFN 
system is an all-or-nothing phenomenon: cells switch on 
signalling or gene expression or not, and the switching cells do 
so at widely variable time points. Nevertheless, we find that a 
reliable antiviral response is achieved through powerful 
paracrine propagation of the signal. Thus, our results show 
that the functional dynamics of the IFN response must be 
understood in terms of the collective spatio-temporal 
dynamics of stochastically reacting single cells. 

Results 

Cell-to-cell heterogeneity in IFN induction 

In order to monitor authentic IFN- (3 expression in real time, we 
stably transfected murine fibroblasts with a BAC-encoded 
reporter expressing TurboGFP under the control of the Ifnb 
promoter (IFN-(3-tGFP). A representative cell clone was 
selected showing stable expression of the reporter. These cells 
were infected with Newcastle Disease Virus (NDV), which 
replicates and induces IFN in the cells via the double-stranded 
RNA sensor RIG-I (Kato etal, 2005), without viral interference 
with this pathway (Childs et al, 2007) . As the newly generated 
viral particles cannot re-infect the mouse cells (Rott, 1979), 
this system allows us to study in a controlled manner the IFN 
induction elicited by the primary infection. To quantitatively 
determine the kinetics and dose response of IFN-(3-tGFP 
expression, reporter cells were infected with NDV and 
subjected to flow cytometry (Figure 1A). We observed a rise 
in IFN-(3-tGFP-positive cells that faithfully reflected the 
accumulation of type-I IFNs in the supernatant (Figure IB). 
The fraction of IFN producers increased nearly linearly over a 
broad range of NDV dose, whereas the average expression 
level already reached ~ 70% of its maximal value at very low 
NDV dose (Figure 1C). The frequency of IFN- (3 expression was 
very similar for different clones, showing the highly repro- 
ducible behaviour of the BAC-based IFN-|3-tGFP reporter 
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(Supplementary Figure S2). Thus, the production of IFN-(3 in 
response to viral infection is controlled by the fraction of 
responding cells. 

To examine whether IFN expression correlates with viral 
replication, we jointly measured the viral protein, hemagglu- 
tinin-neuraminidase (HN), and IFN-(3-tGFP by flow cytometry. 
In agreement with previous findings (Kumagai et al, 2009; 
Rehwinkel et al, 2010), only cells with replicating virus 
expressed IFN-(3-tGFP. However, a large proportion of cells 
with replicating virus did not activate the Ifnb promoter. 
Strikingly, we observed no correlation between the extent of 
replication and the fraction of IFN-(3-tGFP-expressing cells 
(Figure 2A and B; Supplementary Figure S3). These observa- 
tions suggest that the presence of replicating virus in a cell is 
necessary but not sufficient to induce IFN-(3. 

As an alternative explanation of heterogeneous IFN expres- 
sion, it has been suggested that defective viruses are primarily 
responsible for inducing IFN during parainfluenza virus type 5 
(PIV5) infection while during normal replication of PIV5 
effective pathogen-associated molecular patterns (PAMPs) are 
not produced or exposed (Killip et al, 2011). Therefore, we 
examined the intracellular RNA of non-responding virus- 
infected cells for its ability to induce IFN in naive cells 
(Figure 2C; Supplementary Figure S4). RNA from non- 
responding cells induced IFN to a comparable extent as RNA 
from responding cells. This finding shows that the viral RNA in 
IFN-producing and non-producing cells is equally capable of 
inducing IFN and thus the heterogeneity of IFN- [3 induction 
occurs despite the presence of IFN-inducing viral RNA. 

Time-lapse microscopic data revealed that the onset of IFN- (3 
expression after infection of cells with NDV varied strongly 
among the IFN- (3 producers, with cell-to-cell differences as 
large as 20 h (Figure 3 A; Supplementary Figure Sll; 
Supplementary Movie SI). While the increase in viral load 
results in earlier onset, the relative temporal variability of IFN- 
P-tGFP expression differs only slightly. To dissect whether this 
variability is due to variable time points of infection or reflects 
cell-intrinsic properties, we bypassed viral infection by 
liposome-transfecting the cells with the dsRNA analogue poly 
I:C (liposome-free delivery of poly I:C did not lead to IFN-(3 
expression) . Among the IFN producers and independent of the 
poly I:C concentration, the onset time of IFN-(3 expression 
varied strongly and was quantitatively comparable to the case 
of viral infection, as seen by the same order of magnitude of the 
coefficient of variation (CV) (Figure 3A). Furthermore, 
unchanged temporal variability was obtained when synchro- 
nizing viral entry by a temperature shift (Supplementary 
Figure S5). Our data show that the IFN-(3 gene is induced with 
widely varying time delays in the producing cells that are not 
due to variable infection times. In summary, we conclude that 
the cell-to-cell heterogeneity in IFN-J3 expression is predomi- 
nantly of cellular origin. 



Both virus-induced signal transduction and Ifnb 
gene expression are sources of heterogeneity 

To analyse mechanistically how cell-to-cell heterogeneity in 
IFN induction arises, we monitored the activation of the key 
transcription factors NF-kB and IRF-7 in dual reporter cells, 
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Figure 1 Quantitative and temporal heterogeneity of IFN-p induction. A BAC-based reporter construct in which the IFN-p gene is replaced by TurboGFP was 
integrated into murine NIH3T3 fibroblasts. A cell clone with a stable integration of the BAC and representative response towards NDV infection was used (error bars 
represent ± s.d. of triplicates). (A) Induction of IFN-p-tGFP expression upon NDV infection. IFN-p reporter cells were infected with NDV for 1 h. Expression of tGFP was 
determined 24 h post-infection by flow cytometry. Representative dot plots are shown for 10, 20, und 40 HAU/ml NDV. (B) IFN-p reporter reflects endogenous IFN 
production. IFN-p-tGFP expression frequencies after infection with 40 HAU/ml NDV were detected at various time points by flow cytometry. Frequencies were plotted 
against time points post-infection (black circles) and compared with titres of type-l IFN in the supernatant (grey rhombs). (C) IFN-p expression frequency increases with 
viral titre. Reporter cells infected with increasing concentrations of NDV (HAU/ml) were subjected to flow cytometry 24 h post-infection. Frequency of IFN-p-tGFP 
expression (circles) following infection with 1, 2, 5, 10, 20, 40, 80, and 100 HAU/ml and the geometric mean of their fluorescence intensity (triangles) are presented. 
Source data is available for this figure in the Supplementary Information. 



representing the IKKa/P/y-NF-KB and IKKs/TBK-l-IRF-3/ 
IRF-7 pathways downstream of the dsRNA sensor RIG-I. 
Fluorescent protein- tagged NF-kB and IRF-7 (p65-YFP and 
IRF-7-CFP) were localized predominantly in the cytoplasm in 
uninfected cells and accumulated in the nucleus after NDV 
infection (Figure 3B). Nuclear translocation of the two 
transcription factors occurred at the same time in a given cell, 
but the joint translocation time varied widely between cells 
(Figure 3C). This finding was corroborated by antibody 
staining of endogenous IRF-3 and NF-kB (p65), showing that 
the activation of both factors upon NDV infection or poly I:C 
stimulation correlated in individual cells (Supplementary 
Figure S6). Thus, strong cell-to-cell variability arises in the 
shared upstream activation pathway of NF-kB and IRF-7. 

To relate IFN induction to the activation of the latent 
transcription factors we used dual reporter cells expressing the 
IRF-7-CFP fusion protein and the transcriptional reporter IFN- 
(3- tGFP (Figure 3D; Supplementary Figure S7). We found that 
the majority of cells translocating IRF-7 to the nucleus also 
activated the Ifnb promoter (91% at 80 HAU/ml NDV) 
(Figure 3E). Only few cells (9%) did not express IFN-p-tGFP 
upon IRF-7 nuclear accumulation, showing the same distribu- 
tion of IRF-7 translocation times as the IFN-(3-expressing cells 
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(Figure 3E, bottom part). No IFN-p-tGFP induction was 
observed without preceding IRF-7 nuclear accumulation. 
Thus, IRF-7 nuclear translocation is strictly required for IFN- 
(3 expression. The tight coupling of the two events also shows 
the validity of the reporters used. In agreement with the results 
in Figure 3C, there was strong temporal heterogeneity in IRF-7 
translocation between the cells (signalling delay from viral 
infection to IRF-7 translocation T sig = 11.7 ±4.0 h). The inter- 
val between IRF-7 translocation and the onset of IFN-p-tGFP 
gene expression also varied considerably (gene expression 
delay, T gen = 3.4 ± 1.5 h; note that the CVs are similar, a sig / 
r sig = 0.34, G gen /r gen = 0.44). These quantitative data show 
that both cytoplasmic signalling from viral entry to the 
activation of latent transcription factors and induction of 
IFN-p expression cause strong heterogeneity in IFN-p 
production. 



Sister-cell analysis reveals cell-intrinsic 
stochasticity 

Cell-to-cell variability can arise from the intrinsic noise of 
biochemical reactions and from extrinsic factors, such as 
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Figure 2 Viral replication is necessary but not sufficient to induce IFN-p expression. (A) Fractional IFN-p expression among productively infected cells. Reporter cells 
were infected with 40 HAU/ml NDV for 1 h. IFN-p-tGFP reporter expression and intracellular NDV HN protein was measured by flow cytometry at indicated time post- 
infection. Dot plots show IFN-p-tGFP expression among productively infected (NDV HN + ) cells at indicated time post-infection. (B) Separate kinetics of viral replication 
and IFN-p expression. Frequency of IFN-p-tGFP (black circles) and NDV HN expression (grey squares) over time. (C) Unresponsiveness is not caused by the absence 
of inducing viral RNA. NDV-infected (80 HAU/ml) IFN-p-tGFP reporter cells were separated into GFP + and GFP~ fractions. Total RNA was isolated and transfected 
into naive IFN-p-tGFP reporter cells (lower graphs). RNA from non-infected cells served as a control (upper graph). The frequency of IFN-p-tGFP-expressing cells 20 h 
after transfection is presented. Source data is available for this figure in the Supplementary Information. 



differences in cell-cycle stage or cellular environment (Elowitz 
etal 2002; Maheshri and O'Shea, 2007; Paixao et al, 2007; Raj 
and van Oudenaarden, 2008; Snijder and Pelkmans, 2011). To 
minimize extrinsic cell-to-cell differences, we analysed sister 
cells after division (Spencer etal, 2009) . Cells that divided after 
the 1 h period of infection were followed. IRF-7 activation in 
sister cells occurred mainly asynchronously, differing by > 2 h 
in -50% of cell pairs (Figure 4A). The time between IRF-7- 
CFP activation and IFN-p-tGFP expression correlated even less 
in sister cells (Figure 4B), consistent with the previously 
described stochastic transcription of the Ifnb gene (Apostolou 
and Thanos, 2008) . The coefficient of determination r 2 was 0.6 
for r sig and 0.34 for T gen , indicating that 40% of the variability 
in signalling and 66% of the variability in Ifnb gene expression 
are uncorrected between sister cells and thus provide an 
estimate for the cell-intrinsic stochasticity. To examine 
whether different viral replication kinetics in sister-cell pairs 
are a source of variability, we subjected cells to poly I:C 
stimulation. Also, under these conditions cells displayed 
largely uncorrelated IRF-7 signalling and IFN-p gene expres- 
sion (Figure 4E and F), with similar coefficients of determina- 
tion (r 2 = 0.54 for T sig and 0.11 for T gen ) as with viral infection. 
Taken together, these findings show that cell-intrinsic 
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stochasticity is a strong source of cell-to-cell heterogeneity in 
IFN-(3 expression. 

We checked whether the observed differences between 
sister cells relate to their time of division. For this purpose, we 
plotted the time differences between sister cells (AT sig and 
AT gen ) versus the time elapsed since cell division. The very 
weak correlations argue against strong control of heterogene- 
ity by the cell cycle (Figure 4C and D for NDV infection and 
Figure 4G and H for poly I:C stimulation). 

Taken together, the sister-cell analysis indicates a role for 
both 'extrinsic' variability between cells and cell-intrinsic 
stochasticity. The cell-intrinsic component is strong, account- 
ing for approximately half of the variability in the kinetics of 
antiviral signalling and IFN-p induction in individual cells. 
This intrinsic stochasticity provides a rationale for the lack of 
correlation of IFN-p expression with the extent of viral 
replication (cf . Figure 2A) . 

Antiviral protection is an IFN concentration- 
dependent switch in individual cells 

The intrinsic stochasticity indicates that the responsiveness 
of IFN-p induction towards virus is not maximized (i.e., 
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Figure 3 Temporal variability in cellular IFN-p induction. (A) IFN-p expression onset in single cells. Variability of response timing is virus-independent. IFN-p-tGFP 
reporter cells infected for 1 h with indicated concentrations of NDV or transfected with poly l:C at given concentrations were subjected to time-lapse microscopy (15 min 
picture intervals). Distribution of tGFP expression onset overtime (scatter plot, n = 456 (NDV), n= 140 (poly l:C)) and CVs are shown. (B) Synchronous activation of 
NF-kB and IRF-7. NIH3T3 cell clone stably expressing the fusion proteins IRF-7-CFP and NF-KB/p65-YFP were infected with 80 HAU/ml NDV for 1 h and subjected to 
time-lapse microscopy. Fluorescence pictures for CFP and YFP were taken every 20 min. Subcellular localization of IRF-7-CFP (left column) and p65-YFP (right 
column) at indicated time after infection. The diagram shows relative nuclear fluorescence for IRF-7-CFP and p65-YFP from sister cells. (C) Synchronicity is 
independent of response time. IRF-7-CFP and p65-YFP initial nuclear translocation were determined in individual cells and plotted against each other (n = 65). 
Coloured dots indicate the frequency of data points. (D) Expression delay of an individual cell. NIH3T3 cell clone stably expressing IRF-7-CFP together with IFN-p-tGFP 
were infected with 80 HAU/ml NDV. Fluorescence pictures for CFP and GFP were taken every 20 min. Subcellular localization of IRF-7-CFP (left column) and IFN-p- 
tGFP expression (right column) at indicated time after infection. Graphs show relative IRF-7-CFP nuclear fluorescence and tGFP intensity. 7" sig : time interval between 
infection and IRF-7-CFP nuclear translocation. r gen : time interval between IRF-7-CFP nuclear translocation and onset of IFN-p-tGFP gene expression. (E) Response 
variation at distinct stages of IFN induction. The starting times for IRF-7-CFP nuclear translocation were plotted against the times of IFN-p-tGFP expression for 
individual cells (n = 315). Source data is available for this figure in the Supplementary Information. 



many cells with replicating virus do not express IFN-p or do so 
only very late after infection). Therefore, we asked how IFN-p 
production translates into antiviral protection. We chose 



IRF-7 as a prototypical ISG (Honda et al, 2005) and measured 
its expression upon IFN-p stimulation in cells stably trans- 
fected with a BAC-encoded IRF-7-mCherry fusion gene 
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Figure 4 Temporal variability of signalling events in sister cells reveals stochasticity. NIH3T3 cell clone stably expressing IRF-7-CFP or IRF-7-TagRFP together with 
IFN-p-tGFP were infected with 80 HAU/ml NDV for 1 h (A-D, n = 38 sister-cell pairs) or transfected with poly l:C (5 u.g/ml) (E-H, n = 36 sister-cell pairs) and subjected 
to time-lapse microscopy (20 min interval). Sister-pair analysis was carried out for IRF-7 nuclear translocation and IFN-p expression onset. Coloured dots indicate the 
frequency of data points. (A, E) Time of IRF-7 nuclear translocation onset (T sig ) in sister-cell pairs. (B, F) Time intervals between IRF-7 nuclear translocation and IFN-p- 
tGFP expression onset (Tg e n) among sister cells. (C, G) Time elapsed from cell division to IRF-7 nuclear translocation (T si g - T di v) of sister cells. (D, H) Time elapsed from 
cell division to IFN-p-tGFP expression (T sig + T gen - r div ). Source data is available for this figure in the Supplementary Information. 



(Supplementary Movie S2). Flow-cytometric analysis 
showed a digital pattern of IRF-7-mCherry levels, with 
distinct expressing and non-expressing subpopulations, 
where the expressing subpopulation increased with 
extracellular IFN-p concentration (Figure 5A). The binary 
IRF-7 response was consistently found for several IRF-7- 
mCherry clones (Supplementary Figure S8); it was not 
caused by competition of cells for IFN-p, as IFN-p was still 
detectable in the supernatant for >30h (Supplementary 
Figure S9). 

To correlate IRF-7-mCherry with the expression of 
endogenous ISGs, we separated IRF-7-mCherry expressing 
and non-expressing cells by fluorescence-activated cell 
sorting (FACS) and subjected isolated RNA from these cells 
to qPCR analysis. All tested ISGs showed enhanced mRNA 
levels in IRF-7-mCherry expressing cells (Figure 5B). These 
findings suggest that a distinct subpopulation of cells that is 
identified by IRF-7-mCherry expression coordinately 
expresses an antiviral gene programme. To test this, we 
stimulated IRF-7-mCherry reporter cells with IFN-p for 
8 or 24 h to allow induction of ISGs and then infected the 
cells with NDV. Virus replicated in IRF-7-mCherry non- 
expressing cells but not in the expressing cells (Figure 5C). 
This result shows that, as predicted, only the IRF-7-mCherry 
positive subpopulation of cells has been protected by IFN. 
Taken together, these data demonstrate that antiviral protec- 
tion is a stochastic, IFN-p concentration-dependent, switch in 
individual cells. 



Mathematical modelling predicts paracrine 
propagation of the IFN response 

Our data show a remarkable heterogeneity of type-I IFN 
expression and subsequent induction of an antiviral state. To 
understand the resulting dynamics of the IFN system in an 
infected cell population, we developed a multi-scale mathe- 
matical model that describes viral replication and the 
expression of IFN and ISGs in individual cells together with 
intercellular communication through secreted IFN (Figure 6A; 
Supplementary Text SI). According to the sister-cell data, 
signalling and gene expression contribute strong cell-intrinsic 
stochasticity (cf . Figure 4) . For this reason, we attribute cell-to- 
cell variability in the model to intrinsic stochasticity in both 
kinds of processes. 

Two sources of heterogeneous IFN induction are the (i) 
activation of the latent NF-kB and IRF transcription factors 
downstream of viral sensing by RIG-I and (ii) induction of the 
Ifnb gene by nuclear NF-kB and IRFs. Importantly, nuclear 
transcription factor concentration and IFN-p protein increase 
quantitatively very rapidly once nuclear translocation and 
gene induction, respectively, have been triggered. By contrast, 
the waiting times for these events in an individual cell are 
much longer, so that IFN-p induction is switch-like in time (cf. 
sample trajectories in Figure 3B and D) . Therefore, we describe 
the dynamics of these events as stochastic transitions 
u;^ ± andu;7 ± between discrete cell states (virus-infected, NF- 
KB/IRF-activated, IFN-p expressing along the horizontal axis 
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Figure 5 Bimodal antiviral response towards IFN. A BAC-based reporter construct containing mCherry fused to the C-terminus of the chromosomal IRF-7 gene was 
integrated into NIH3T3 cells. Experiments were performed with a cell clone exhibiting a stable integration of the BAC and a representative response towards IFN. 
(A) Binary dose- and time-dependent IRF-7-mCherry expression. IRF-7-mCherry reporter cells were stimulated with the indicated concentrations of IFN-p, mCherry 
expression was determined by flow cytometry. (B, C) Bimodality of IRF-7 expression is reflected in ISG transcription and antiviral protection. (B) Reporter cells were 
treated with 500 U/ml IFN-p for 16 h and subjected to FACS. IRF-7-mCherry positive (2) and negative (1) populations indicated by the shaded areas were separated. 
RNA was prepared from both populations and analysed by qRT-PCR for expression of the indicated ISGs and Rps9 as a control. The results were normalized to (3-actin 
mRNA and are shown as fold increase of untreated reporter cells. (C) Reporter cells treated with 500 U/ml IFN-p for 8 and 24 h were infected with 80 HAU/ml NDV. In all, 
20 h after infection, cells positive (blue) or negative (red) for NDV HN (left column, intracellular antibody staining) were analysed for IFN-stimulated IRF-7-mCherry 
expression (right column) by flow cytometry. 
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in Figure 6A) (Mariani et al, 2010) . The intracellular replication 
of virus particles is modelled as a birth-death process. 
Secreted IFN diffuses in the extracellular space and induces 
ISGs in an autocrine or paracrine manner. Given the high 
diffusion coefficient of IFN (Kreuz and Levy, 1965), diffusion 
on the relevant length scale of several cell diameters takes 
place within minutes, compared with the hour-timescale of 
IFN and ISG expression, so that we assume uniform distribu- 
tion of secreted IFN. To account for the switch-like induction of 
ISGs, such as IRF-7 (cf. Figure 5), we modelled the IFN- 
dependent activation rate of STAT1/2 nuclear translocation, 
Wstat^ (Lillemeier et al, 2001; Vinkemeier, 2004), followed by 
the rate of ISG induction, resulting in a bimodal 

distribution of ISG-expressing and non-expressing cells. Virus 
ceases to replicate in cells expressing ISGs. 

To determine the model parameters from the experimental 
data, we focused on our extensive data for high-dose 
infections, 40 and 80 HAU/ml NDV (Figure 6; Supplementary 
Text SI; Supplementary Figure S10). We first computed the 
distributions for the NF-kB/IRF-7 nuclear translocation times 
and the IFN-(3-tGFP switching-on times. Parameter optimiza- 
tion using as objective function the squared sum of residuals 
between measured and computed (binned) distributions 
yielded a good fit (Figure 6B and C and green trace in 
Figure 6D). This fit required that (i) induction of RIG-I- 
mediated signalling by virus is cooperative (Supplementary 
Figure S10B) and (ii) Ifnb gene induction is modelled as a 
multi-step process (Supplementary Figure S10C). Both impli- 
cations of the quantitative modelling of the data — cooperative 
RIG-I signalling and multi-step assembly of the IFN- (3 
enhanceosome — are consistent with experimental findings 
(Ford and Thanos, 2010; Onoguchi et al, 2010). 

Next, we parameterized IRF-7/ISG induction by fitting the 
dose response of IRF-7-expressing cells versus external IFN-(3 
(Supplementary Figure S10D) and determining the IFN- (3 
secretion rate per cell to match the observed time course of 
extracellular IFN- (3 (Figure 6D, blue trace) to the fraction of 
IFN-(3 producers (Figure 6D, green trace). With these 
parameters (Supplementary Table SI), the simulated kinetics 
of IRF-7-expressing cells agreed with the data (Figure 6D, red 
trace; note that IRF-7 induction takes place during prolonged 
elevation of extracellular IFN-(3). Importantly, the stochastic 
switching in individual cells translates into predictable 
dynamics of IFN-secreting and protected cell fractions at the 
population level. 

The comparison between model and data in Figure 6D 
suggests that the IFN- [3 produced by a single cell can induce 
ISGs in several cells. To examine this further, we simulated the 



model for smaller infection doses where IFN- (3-pro during cells 
are sparser. Measuring ISG expression, we found good 
agreement between model and experiment over a range of 
small doses of infection (Figure 7A). Importantly, the model 
predicted strong paracrine amplification of antiviral protection 
through IFN, with up to 40 times as many ISG-expressing cells 
as IFN-secreting cells. This prediction proved remarkably 
accurate when tested experimentally (Figure 7B and C). Thus, 
at the level of the cell population, paracrine amplification 
transforms stochastic single-cell responses into efficient and 
predictable antiviral protection. Few IFN-producing cells 
suffice as sentinels of viral infection to protect a large number 
of cells surrounding the infection site. 



Discussion 

In this paper, we have dissected the dynamics of the type-I IFN 
system in an experimental setting that allowed reliable 
quantitation in living cells. In addition to heterogeneous IFN 
induction, we describe for the first time bimodal, all-or- 
nothing expression of ISGs that governs resistance to viral 
infection. Thus, heterogeneity of single-cell decisions is a 
pervasive feature of the IFN system. Imaging signal transduc- 
tion and gene-expression dynamics, we show that these 
processes contribute quantitatively to heterogeneous expres- 
sion of the Ifnb gene, and that cell-intrinsic stochasticity 
prevails in both antiviral signalling and IFN- (3 expression. 
Importantly, however, we find that paracrine propagation of 
the IFN signal can mediate reliable antiviral protection at the 
cell-population level. 

The origin of heterogeneous IFN induction has been 
controversial; both cell-intrinsic stochastic events in gene 
expression (Hu et al, 2007; Apostolou and Thanos, 2008) and 
virus-related mechanisms (Chen et al, 2010; Killip et al, 2011) 
have been implicated. Our live-cell imaging data obtained with 
both virus and viral-surrogate (poly I:C) stimuli demonstrate 
that cell-intrinsic mechanisms generate large-scale heteroge- 
neity in the decision to express IFN- (3 and in the timing of 
expression. This establishes intrinsic host cell variability as a 
decisive factor shaping the antiviral response at the single-cell 
level while not excluding additional virus effects in other 
systems (e.g., interference with IFN-p induction, which is a 
hallmark of many pathogenic viruses) . 

Cell-intrinsic stochasticity occurs at multiple levels of 
regulation to generate heterogeneity in IFN- (3 expression. 
Quantifying the relative contributions of gene expression and 
signal transduction activating IRFs and NF-kB, we found that 



Figure 6 A multi-scale mathematical model of IFN induction and response. (A) Scheme of the state transitions of an individual cell. The model describes a population 
of individual cells with virus replication, induction of IFN genes through RIG-I signalling to NF-kB and IRFs, induction of ISGs (including IRF-7) by IFN, and cell 
communication via secreted IFN. The possible state transitions of an individual cell with the propensity functions ware shown. The colour code indicates nuclear 
translocation of IRFs/NF-kB (cyan nucleus), induction of IFN (green cytoplasm), and induction of IRF-7 (red cytoplasm), as imaged experimentally. Induction of IFN-p 
(dashed arrow W\ + ) is explicitly modelled as a multi-step process to match the available experimental data (see (C) below); all other transitions are modelled as single 
steps (solid arrows). (B) Model simulation of NF-kB/IRF nuclear translocation times versus experimental data. The distribution computed with the model (black line) 
matches the data (histogram) obtained for high infection dose (80 HAU/ml). (C) Model simulation of IFN-p-tGFP onset times versus experimental data. The distribution 
computed with the model (black line) matches the data (histogram) obtained for high infection dose (80 HAU/ml). (D) Comparison between model and experimental data 
for high-dose NDV infection. The model reproduces the observed dynamics of viral load (as measured by HN expression), IFN-p-tGFP induction, extracellular IFN titre 
and IRF-7-mCherry expression (solid lines, model; dots, experimental data; infection dose 40 HAU/ml). Note that the smooth model curves in (B, C) are mean values 
obtained by simulating 10 4 cells. Source data is available for this figure in the Supplementary Information. 
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both processes have strong stochastic components. The 
decision whether to express the Ifrib gene was made primarily 
at the level of signal transduction, indicating a limitation in 



RIG-I-induced signalling. Stochasticity in gene expression 
contributed to temporal variability. These data show that 
heterogeneous IFN- (3 expression cannot be attributed solely to 
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Figure 7 Paracrine amplification of the IFN response is predicted by the model and verified experimentally. (A) Comparison between model and experiment for low 
infection doses. For low viral doses yielding sparse infection (1 , 2, and 5 HAU/ml), the model predicts IFN response kinetics of the correct magnitude (solid lines, model; 
dots, experimental data). (B) Paracrine propagation of the signal. The rarity of IFN-p-producing cells at low infection doses reveals strong paracrine propagation of the 
IFN response because up to ~40 times as many cells respond with IRF-7-mCherry expression (IFN responders). The agreement between model (squares) and 
experimental data (black dots) is remarkable, given that the model was calibrated only for a single condition (40 HAU/ml, marked with blue square) while the red squares 
represent predictions that were subsequently tested. (C) Co-culture of IFN-p-tGFP (green) and IRF-7-mCherry (red) cells illustrates paracrine communication. IFN-p- 
tGFP reporter cells were infected with 40 HAU/ml NDV for 1 h. Then IRF-7-mCherry cells were added at same density and cells were subjected to time-lapse 
microscopy. Merged fluorescent pictures for IFN-p-tGFP and IRF-7-mCherry at selected time points post-infection are shown. Source data is available for this figure in 
the Supplementary Information. 



stochastic gene regulation, which has previously been the 
main suspect for IFN- (3 and other cytokines (Hu-Li et al, 2001; 
Guo et al, 2004; Apostolou and Thanos, 2008; Raj and van 
Oudenaarden, 2008; Mariani et al, 2010). Rather, both virus- 
induced signalling and gene regulation could have different 
control on cell-to-cell variability in different cell types, 
depending, for example, on the expression levels of pathway 
components. 

Our model prediction that viral induction of signalling 
exhibits cooperativity might be linked to the clustering of the 
adaptor IPS-1 downstream of RIG-I (Onoguchi et al, 2010). 
Indeed, the formation of large macromolecular complexes in 
signalling transduction as well as gene expression (Ford and 
Thanos, 2010) are likely to be slow stochastic processes 
(Luijsterburg et al, 2010) that could determine both frequency 
and strong temporal heterogeneity of IFN- (3 expression. 
Distributed expression levels of signalling molecules, tran- 
scription factors etc. may further add to this cell-intrinsic 
stochasticity in creating cell-to-cell heterogeneity (Feinerman 
et al, 2008) . A recent study has correlated variability of both 
viral and cellular factors with stochastic IFN-p expression 
(Zhao et al, 2012). 
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Despite the large heterogeneity of single-cell decisions to 
express IFN- (3 and ISGs, the cell-population response is 
remarkably predictable. Our data-driven model indicates that 
this is a robust feature of the system that requires simply an 
appropriate tuning of IFN-p production per cell and IFN 
responsiveness of ISG induction, allowing paracrine propaga- 
tion of the IFN signal. In this way, an individual IFN-secreting 
cell can bring a defined territory surrounding it into an 
antiviral state. 

Thus, heterogeneous and stochastic IFN induction could 
primarily reflect mechanistic limitations of the host cells (as 
has been speculated for other cytokines; Mariani et al, 2010) 
that are readily overcome through paracrine signal propaga- 
tion. Alternatively, heterogeneity may have functional bene- 
fits. In particular, the widely distributed onset times of 
producers could sustain signalling if the IFN production period 
per cell is limited. Moreover, we expect that different 
interference strategies of pathogenic viruses, inhibiting IFN 
versus ISG induction (Versteeg and Garcia-Sastre, 2010), will 
have differential effects on the spatio-temporal dynamics of 
the IFN system which, in turn, might control the viral 'choice' 
between persistent versus acute lytic infection. 
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Our demonstration of paracrine cell-to-cell communication 
as an efficient mechanism for quenching single-cell stochas- 
ticity will likely be relevant also for other systems. It 
underscores the importance of studying single cell behaviour 
within cell-population context for rationalizing cell 
physiology. 

Materials and methods 
Plasmids and BAC constructs 

Cloning of BACs was performed by homologous recombination (Zhang 
etal, 2000) in Escherichia coli strain GS 1783 (Tischer et al, 2010). BACs 
of the RPCI-24 library were obtained from Children's Hospital Oakland 
Research Institute. IFN-(3-tGFP is based on clone RP24-325J11, where 
the Ifnb ORF was replaced by TurboGFP ORF (Evrogen) containing an 
N-terminal triple FLAG epitope and a FRT-site flanked neomycine 
phosphotransferase gene under control of the PGR promoter as 
described earlier (Pulverer et al, 2010). IRF-7-mCherry is based on 
clone RP24-266F8, in which the Irf7 stop codon and 3' adjacent 
~ 300 bp were replaced by the mCherry (Shaner et al, 2004) coding 
sequence and the above-mentioned neomycine selection cassette. 
Cassettes containing TurboGFP or mCherry together with the selection 
marker were amplified by PCR introducing 50 bp target sequence 
homologies (Ifnb-tGFP forward: 5'-GCCTTTGCCTCATCTTGCAGGTA 
GCAGCCGACACCAGCCTGGCTTCCATCATGGACTACAAAGACCACGAC 
GGAG-3'; Ifnb-tGFP reverse: S'-GATGCCACAGTCATTTCCCTAGAT 
GCTTTAGATTTCCCATTATCTCTTACCCTTACAATTTACGCCCTTAAGAT 
CC-3'; Irf7-mCherry forward: 5'-CAACAGTCTCTACGAAGACATCGAA 
CACTTCCTCATGGACCTGGGTCAGTGGCCTTCAGGTGGAGGCGGTATG 
GTGAGCAAGGGCGAGGAGGATA-3'; Irf7-mCherry reverse: 5'-CCCC 
CAGTCCCTACACACCAGGCCTCCTTTGCAGCACTCACTTTCATAATCAA 
GACTCAAGGGCATCGGTCGACGGATCC-3'). PCR products were trans- 
ferred into E. coli by electroporation. Correctness of recombination and 
integrity of the BAC molecules was checked by PCR, restriction 
analysis and sequencing of the integration site. 

p65-YFP was created by inserting amplified murine p65 cDNA into 
pMBC-1 (Dirks et al, 1994) containing linker and YFP encoding 
sequences. pIRF-7-CFP and pIRF7-tagRFP were created by introdu- 
cing Irf7 cDNA (C57BL/6) (amplified by forward: 5'-CGAATTCCACCA 
TGGCTGAAGTGAGG-3'; reverse: 5'-CGCCGGCCCGCCTCCACCTGAAG 
GCCACTGACCCAGGTCC-3') into pMBC-1 vector containing linker and 
sequences for CFP and tagRFP (Evrogen), respectively, via EcoRI and 
Nael restriction sites. 



Cell culture and DNA transfection 

Murine fibroblastoid NIH 3T3 cells were grown under standard tissue 
culture conditions in DMEM medium (Sigma) supplemented with 10 % 
fetal calf serum, glutamine, penicillin and streptomycin. Transfections 
of BAC or plasmid DNA were performed with Metafectene (Biontex) 
according to the manufacturer's instructions. G418- or puromycin- 
selected representative clonal cells showing stable expression of the 
reporter construct and strong signal to background ratio of the 
fluorescent marker were used. Transfections of BAC or plasmid DNA 
were performed with Metafectene (Biontex) . 



IFN stimulation 

Infections with Newcastle Disease Virus LaSota (Lohmann Tierzucht, 
Cuxhaven, Germany) were performed after three times washing with 
serum-free medium. After 1 h, residual virus was removed by washing 
three times with serum-containing medium. Low molecular weight 
poly I:C (InvivoGen) was transfected via Lipofectamine2000 (Invitro- 
gen) according to the manufacturer's instructions. Total RNA of 
fractionated cell populations was isolated using RNeasy (QIAGEN) 
according to the manufacturer's instruction and transfected using 
Lipofectamine2000. 



IFN source and detection 

Murine IFN-p was recombinantly produced in stably expressing BHK- 
21 cells and obtained from supernatant as described before (Pulverer 
et al, 2010). IFN activity was measured with a bioassay based on Mx2- 
Luc reporter gene expression (Kugel et al, 2011). 



Time-lapse microscopy 

Time-lapse microscopy used an Olympus 1X81 fluorescence micro- 
scope, with a PRIOR motorized stage and climate chamber. The 
Olympus Cell A M autofocus function was used to compensate Z-drift. 
Cells were seeded into ibidi eight-well u-slides coated with collagen IV. 
Image acquisition was performed at intervals between 15 and 30 min. 



Image and data analysis 

Microscopical picture series were analysed with ImageJ (NIH, 
Bethesda, MD) and built-in plugins as well as MTrackJ (E Meijering). 
Time points of p65 or IRF-7 activation were plotted upon significant 
increase of nuclear signal that eventually led to maximal accumula- 
tion. Time of cell division was counted when sister cells were clearly 
separated. Figures and statistical analysis were done with MATLAB, 
Microsoft Excel 2007 and GraphPad Prism 5. 



Immunocytochemistry and flow cytometry 

NDV-infected cells were stained in suspension following fixation with 
4% paraformaldehyde (10 min) and Triton X-100 (0.1%) permeabili- 
zation. Primary mouse anti-HN14f (Santa Cruz) was used at 1:500. 
After washing with PBS supplemented with saponin (0.1%), 
secondary Cy3- or Cy5-labelled antibody (goat anti-mouse IgG + IgM, 
Dianova) was used. Analysis of living and fixed cells in flow cytometry 
was performed with BD FACscalibur and LSR II (analytical) or 
FACSaria (sorting), using appropriate laser and filter settings for 
TurboGFP, mCherry, and chemical fluorophores. Results were quanti- 
fied with FlowJo 7.6 software. 



Real-time PCR 

mRNA was isolated with RNeasy (QIAGEN), cDNA hybrid strands 
were generated with first-strand cDNA synthesis kit (GE Healthcare). 
Reactions were performed with QIAGEN QuantiTect SYBR green kit in 
a Roche LightCycler 480 using specific intron-spanning primers (Mx2: 
forward 5' -TCACCAGAGTGCAAGTGAGG-3 ' , reverse 5'-CATTCTCCCTC 
TGCCACATT-3'; Rsad2: forward 5' GTCCTGTTTGGTGCCTGAAT-3', 
reverse 5'-GCCACGCTTCAGAAACATCT-3'; Isg20: forward 5'-TAAGCG 
CCTGCTACACAAGA-3', reverse S'-GCAGCTTCTAACCCTGGATG-S'; 
Uspl8: forward S'-AAGGACCAGATCACGGACAC-S', reverse 5'-CATC 
CTCCAGGGTTTTCAGA-3'; Irf9: forward 5 ' - ACCACGGAACCAGAAATC 
AC-3', reverse 5'-GTTGCAGTTGCTGTTGCTGT-3'; Irf7: forward 5'-GAA 
GACCCTGATCCTGGTGA-3', reverse 5'-CCAGGTCCATGAGGAAGT 
GT-3'; ML: forward 5 '-CTCACCAGGAACCAGAGGAA-3 ' , reverse 5'-TG 
AGTGGTGTAACTGCTGTGG-3'; Actb: forward 5'-TGGAATCCTGTGGC 
ATCCATGAAA-3', reverse 5'-TAAAACGCAGCTCAGTAACAGTCCG-3'). 



Mathematical modelling 

We simulated in 10 4 cells the stochastic intracellular dynamics for virus 
replication, IRF and NF-kB activation by virus, induction of IFN, 
activation of STAT1/2 by IFN, and expression of IRF-7 as a 
representative ISG, using the Gillespie's algorithm. Extracellular IFN 
was considered as rapidly diffusing in cell culture and iterated in 
parallel with a deterministic Euler method. Initially, signal transduc- 
tion and gene expression were inactive in all cells, and the cells were 
infected with virus. Parameters were chosen according to experimental 
measurements as explained in Supplementary Information. MATLAB 
codes for all model simulations are available as Supplementary 
Information. 
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Supplementary information 

Supplementary information is available at the Molecular Systems 
Biology website (www.nature.com/msb). 
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